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Abstract 

We describe the light-cone wavefunction in the saturation regime in terms 
of the density of gluons per unit of transverse phase space, the occupation 
number, and in terms of the color charge correlator. The simple McLerran- 
Venugopalan model gives what are claimed to be general results for the phase 
space gluon density, but it does not well describe the general case for the 
charge correlator. We derive the general momentum dependence and ra- 
pidity dependence of the color charge correlator which exhibits strong color 
shielding. A simple physical picture which leads to these general results is 
described. 



1 Introduction 

Our object in this paper is to describe what the light-cone wavefunction looks 
like in the saturation regime, that is when the occupation numbers of gluons 
are large. What we attempt to describe is the density of gluons and the color 
charge correlator as a function of rapidity and transverse momentum. We 
shall carry out this description first for the McLerran-Venugopalan model [ttj 
|2|. [|, then for an unsaturated BFKL[|, Q light-cone wavefunction and finally 
for a more general saturated wavefunction, both for fixed coupling and for 
running coupling. 

lr This research is supported in part by the US Department of Energy. 
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For the gluon density in phase space we find the result that rf a b ^ is 

proportional to JV 'C 1 £n[Q 2 s {Y) /IjJ in all the saturation models. We believe 
this to be a very general result which follows straightforwardly from the 
physical picture we give in Sec. 6. 

The color charge correlator p(£±, Y), defined below in (17), is a constant 
in the McLerran-Venugopalan model reflecting the fact that the elementary 
charges are completely uncorrelated in that model on scales small compared 
to one fermi. We note that the first gluonic corrections to the model, cal- 
culated in Sec. 2 begin to show charge shielding. Charge shielding is clearly 
visible in the BFKL wavefunction and in saturated wavefunctions produced 
from BFKL evolution. In the latter case, both for fixed coupling and for 
running coupling, -j^p(£±,Y) = Nc ~ X £\ in the saturation region, and the 
proportionality to £\ means that charge shielding is as complete as it can be. 
That is when measured on a scale Ax± ~ l/£± all gluons having momenta 
k± ^> £± contribute not at all to the charge correlator. 

In Sec. 6 we present a simple picture which leads to the results quoted 
for Jbd?i an d fo r ippC^-U ^0- This picture is most easily explained by referring 
to Figs. 5 and 6 where the saturation region is shown. In Fig. 5 the saturation 
region is illustrated for fixed coupling in the y and in k\/A 2 plane. Note that 
y increases in the downward direction in the figure. The diagonal line is the 
separation line between the dilute and dense regimes. In Fig. 6 the saturation 
region is shown as a function of x_ and y in the Kovchegov gauge 0. (This 
picture is that first given by Iancu, Leonidov and McLerranJfj where a gauge 
was used which puts all the gluons in saturation at positive values of x_, 
rather than the choice here of negative values of ||.) In terms of these 
figures, gluons contributing to -^p(£±, Y) are located in the shaded region C 
and they come by direct bremsstrahlung emission from the shaded region A. 
In the region A, £± is the saturation momentum of that rapidity and so the 
sources at A do not have strong shielding on the scale £± and thus the gluons 
they emit into region C, having a comparable transverse momentum £±, are 
also not far from random in color. This naturally gives -j^p(£±,Y) oc £\ 
with the 1/a density of sources being cancelled by the a for a gluon emission 
leaving -^p independent of a. Gluons can come into the region C from 
regions other than A but they do not contribute to -^p because they are 
strongly shielded, as explained in Sec. 6. 

Now turn to i?5L , evaluated at Y. Such gluons come generically from 
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anywhere along the saturation line in Fig. 5, the region B being typical. In 
Sec. 6 we show how emissions from the saturation line, between Y(£±) and Y, 
naturally lead to Jf^- h ~ ± • £n[Q 2 (Y)/£ 2 ^ g g §. 

Finally, in Sec. 6 we show how the x_— distribution of saturated gluons, 
as shown in Fig. 6, naturally leads to a Glauber-like picture of the scattering of 
a color dipole on a saturated wavefunction. In passing the saturated hadron 
the dipole encounters, in a sequential manner, gluons located at different 
rapidities and scatters independently off them thus naturally leading to the 



exponential result required by the Balitsky-Kovchegov equation [0, [L^, [13 



2 The McLerran-Venugopalan model [1-3] 
2.1 The valence quark approximation 

In the leading approximation the wavefunction of a high momentum large 
nucleus consists of the valence quarks of the various nucleons of the nucleus 
Lorentz-contracted to a thin disc. So long as transverse distances less than a 
fermi are considered the quarks may be viewed as uncorrelated in color. At 
an impact parameter b from the center of the nucleus the number of valence 
quarks per unit area is 



N(b) = N N C 2VR 2 - b 2 (1) 

where Nq is the nuclear density, the number of nucleons per unit volume in 
the rest system of the nucleus, N c is the number of colors, and R is the 
radius of the nucleus. Throughout this paper we take a simplistic view of a 
nucleus as having a uniform density, No, of nucleons within a sphere of radius 
R. 

Suppose we take the nucleus to be moving along the positive z-axis. Now 
scatter a left-moving quark-antiquark dipole, moving along the negative z- 
axis, on the nucleus. The probability that the dipole pass through the nucleus 
without interacting is equal to 



S 2 (b,x ± )=e- N ®°*^ (2) 
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Figure 1: 



where x± is the transverse separation of the quark- ant iquark pair of the 
dipole and S is the S-matrix for elastic scattering of the dipole with the 
nucleus. a q is the cross-section for scattering of the dipole on a single quark 



which we take to be [114" 



<y q {xi_) = -j^-x]_xG q (x, l/x\) (3) 
where the gluon distribution in a quark is given as 

xG q (x, 1/xl) = ^£n (4) (4) 

in the two-gluon exchange approximation illustrated in Fig.l. In (4) r$ is the 
radius of a nucleon and we always suppose that x\jr\ 1. The S-matrix 
can be written as 



S = e-^Qll± ( 5 ) 
where the saturation momentum measured by quarks is given as 

Ql = C F /N C Q 2 S (6) 
and where || the gluon saturation momentum is 
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Q^j^N^xG^Ql) (7) 
with the gluon distribution of a nucleoli xG given by 

xG((x,Q 2 s ) = N c xG q (x,Q 2 s ). (8) 



Finally, the color charge correlator in the nucleus takes the form|T5 

< p\b + x x )p>(b) > MV = ^-^5 2 {x A _)g 2 C F N{b) (9) 

c 

where we use a normalization such that 



Q 2 

< P l {b + x L )p>(b) > quark = 5 ij 6 2 (x ± )^- (10) 



for the color correlator of a single high momentum quark. 

2.2 Valence quarks along with Weizsacker- Williams glu- 
ons 

The valence quarks of a fast nucleus can emit (virtual) gluons which then 
form the gluon distribution of the nucleus in the semiclassical, Weizsacker- 
Williams, approximation. The well-known result for the distribution of glu- 



cLxGa 

d 2 e ± d 2 b 

or 



f d 2 x L N 2 - 1 2 Q z /4 ilnx . x± 

J 4tt 4 aN c x 2[ } 1 } 



dH ± cPb £ 2 JQ 2 <^1 47rW £(g ' /ln ^' (12) 
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2(n*-i) cPi^cPb * s ^ e occu P a ti° n number of gluons, per unit of rapidity. We 
believe (12), except for overall normalization, is a general result for the phase 
space density of gluons deep in the saturation regime, ^±/Q 2 <C 1, of a high 
energy wavefunction ||, |TC|j . 

Now consider the scattering of a dipole on the nucleus, including the semi- 
classical gluons in the wavefunction of the nucleus. We wish to determine 
how (5) and (9) change in the presence of the Weiszacker- Williams gluons 
in the nucleus. The calculation is most easily done by Lorentz transforming 
the dipole-nucleus system so that the dipole carries a greater momentum and 
the nucleus a lesser momentum. In this way we may view the Weiszacker- 
Williams gluons involved in the interaction as being part of the wavefunction 
of the dipole. The scattering which modifies (5) now can be taken to be that 
of a left-moving quark-antiquark-gluon system on a McLerran-Venugopalan 
nucleus having only valence quarks in its wavefunction. In the large N c limit 
this contribution, 6S, is given by 

aN 

SS(b,x ± )=AY ' 



2tt 

d 2 z±x 2 ± 



(x 1± - z ± ) 2 (z± - x 2 ±f 



S{b,xi± - z±)S(b, z± - x 2 ±) (13) 



and is illustrated in Fig. 2. Deep in the saturation regime, Q 2 x\ ^> 1, the 
dominant contribution to (13) comes when z± = Xl± + X2± . One finds 



where AY is the range of rapidity of the gluons we are considering. Eq.14 
shows that the transmission amplitude actually grows and becomes greater 
for a quark-antiquark-gluon to pass through the nucleus without interaction 
than for a quark- ant iquark to pass through without interaction. In the frame 
in which we have done the calculation the extra gluon tends to shield the 
charges in the dipole. 
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Figure 2: 

In a frame where the Weiszacker- Williams gluons are in the nuclear wave- 
function, and (14) is given by the gluon dependent part of the scattering of 
a quark- ant iquark dipole on the nuclear wavefunction, the shielding is ex- 
pressed by 

< p*(x ± + b)p>{b) >ww= ^ < PK X ± + &V ( 6 ) >mv ■ (15) 

This shielding continues to grow as more and more gluons are created. The 
McLerran-Venugopalan valence quark model seems to be a good starting 
point for calculating the density of gluons in the saturation region, as given 
by (11) and (12), but it is not a good starting point for evaluating the dis- 
tribution of color charge in a large nucleus. In the McLerran-Venugopalan 
model charges are randomly distributed while in the actual wavefunction 
there is an arrangement of the quarks and gluons so as to give strong color 
shielding. Eq.(15), or (14), is the first indication of that shielding. 

3 The BFKL wavefunction 

The BFKL0, || wavefunction in a dilute (non-saturation) regime already 
exhibits the charge shielding which we shall find in the saturation regime 
at extremely high energies. The gluon distribution in the BFKL wavefunc- 
tion of a high-energy hadron is related to the charge correlator of gluons 
at all higher rapidities by 
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d £ = ^w^n (i6) 



Our notation is such that 



p{£ ± , Y)= J d 2 x ± e- u ^ x± < P l ( b + x ±)P i ( b ) > d2b 



(17) 



and where the Y variable now indicates a dependence on rapidity in contrast 
to the McLerran-Venugopalan model. Using 



d 2 £ ± ^aY e ± 

for fixed coupling BFKL evolution one finds 



dxG e ^P-l)Y 1 

c —f^t ir ( 18 ) 



e (a P -l)Y 

p(£±,Y) =4n 3 c£ ± ^_ . (19) 



The corresponding result using (9) is 



Pmv(£±) = 9 2 C F ■ N q (20) 

with N q the total number of quarks in the nucleus. Aside from the Y- 
dependence in (19), which one expects will disappear in the saturation regime, 
the big difference between (19) and (20) is that puv has no dependence 
on £± while the BFKL correlator grows linearly with £±. The lack of an 
£±— dependence in p MV is natural and expresses the fact that the quarks are 
completely uncorrelated in color in our short distance approximation. The 
growth of the BFKL correlator with £±, as given by (19), indicates strong 
color shielding. The linear growth in £±_ in (19) means that source gluons 
having k± <^ £± contribute little to p(£±, Y) since there are two few such 
gluons while source gluons having k± 3> £± contribute little to p(£±, Y) so 
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that such high momentum gluons are perfectly shielded. Consider an "area" 
in momentum space A£\ = \7i£ 2 ± , then 



d/dYp(£ ± ,Y) (2tt) 2 

A£ 2 ± dxG/d 2 £ ± -— {ap ~ 1} (21) 

being independent of Y and £± indicates that source gluons having k± £ ± 
contribute essentially independently and incoherently to the charge correla- 
tor. 

Thus the picture of the distribution of color charges in a BFKL wave- 
function is quite simple. Gluons having momentum £± appear as completely 
shielded by other gluons having a similar momentum when viewed on a scale 
Ax± > l/£± while they appear as completely random charges when viewed 
on a scale Ax± < l/£± with the transition between the two pictures occuring 
at Ax± ~ l/£±. We shall see that this picture is preserved in the saturation 
region where BFKL evolution no longer applies. 



4 The light- cone wave function in the satura- 
tion regime 

To deal with the saturation regime we use the Kovchegov equation [12]. (The 



Balitsky equation [11] would give the same results.) We take the equation in 



the form, with x± = Xi± — x 2 ± 



dS(x±,b,Y) aCp C d 2 z±x 2 L 

dY 7T 2 J (ar 1JL - z x ) 2 (z± - x 2 ±) 2 

■ [S(x 1± - z ± , b, Y)S(z ± - x 2± , b, Y) - 3{x ± , b, Y)} , (22) 

as illustrated in Fig. 3. We suppose the vertices in the upper parts of the 
graphs are taken in the large N c limit, though of course we cannot assume 
a large iV c limit for all the graphs as unitarity effects would be lost. We 
write Cf rather than N c /2 in the prefactor in (22) because that is exact for 
the virtual correction which will dominate our result. The factorized form of 
the two S's on the righthand side of (22), which was originally motivated by 
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Figure 3: 



Kovchegov by choosing scattering on a large nucleus, cannot be expected to 
be exact, but this will not be important for the argument which follows. 



If Q 2 {Y)x\ ^> 1 the quadratic term in (22) is small when (x\± — z±) 2 Q 2 
and (x2± — z±) 2 Q 2 are both greater than 1. Thus the dominant contribution 
to the righthand side of (22) comes from the second term in the regions 



1/Ql < (x 1± - z ± ) 2 < x\ 



(23) 



and 



l/Ql < (x 2 ± - z ± ) 2 < x\ 



(24) 



giving 



dS(x±,b,Y) 
dY 



2aC F 



TV 



Y). 



(25) 



Using 



HQ 2 S (Y)/Q 2 S (Y )} 



(Y - Y ) 



(26) 



c 



with 0, |T§, |T7| 
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Figure 4: 



C= MM (27) 

1 — An 

where x(^) is the usual eigenvalue function appearing in BFKL evolution, 



and An is determined by x = — vA - - Thus|]13 



S(x ± , b, Y) = exp | ln 2 [Q2( & , y )^] J 5 6j y , 



(28) 



where = l/a^. 

Moving again to a frame where the gluon at z±_ in (22) and in Fig. 3 is in 
the wavefunction of the hadron on which the dipole scatters, it is natural to 
identify the exponent in (28) as the scattering, due to two gluon exchange, 
of the dipole on the saturated hadron. This gives 



cN c ™ ' ' - 2N C J 1/X 2 r-'l-' Itt 

as illustrated in Fig. 4 while Eq.29 gives 



p(4±,6,y) = — £*nft£/£], (30) 

COL 
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where c is given in (27). Now (16) is no longer exactly valid, however, as we 
shall see in Sec. 6, it does continue to have approximate validity. Using (16) 
with an undetermined constant k gives 



dxG . N?-l 1 - A 



dw ± h ikdi i^^^- (31) 

This result is exactly the same as that found in Ref.9 and except for the 
way the constants are written that of Ref.10 and of the original McLerran- 
Venugopalan model, (12). We believe it to be a very general result. 



5 Saturation with a running coupling 

In this section we shall generalize the results of Sec. 4 to the case of QCD 
with a running coupling. We begin with the generalization of (25) to the 
running coupling case which now reads 



dS(x±,b,Y) 
dY 



2C F 
Tib 



in 



lnQ 2 s /K 2 ' 
—inx 2 , A 2 



S(x ± ,b,Y). 



(32) 



The analogue of (26) isJTBI, |TT 



in 



Q 2 s(Y,b) 
A 2 



2Nc 
nb 



cY +o(r 1/6 ) 



(33) 



where c is exactly as given by (27). Combining (33) with (32) yields 



S(x±, b, Y) = exp 



\ -C F 
\2cN c 



^ 2/A2l „, ln^/A 2 A ,^, 2l „ Ql 



2 In 2 Qj /A 2 In 



— In x 2 , A 2 



— \nQ s x ± In 



x 2 A 4 



■S (x±,b,Y(x±)) 



(34) 



which now replaces (28). 

Using the exponent in (34) to identify p, exactly as was done in (29), 
gives 
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p(e±, b,Y) = ^<->£l HQl/ll] (35) 
c a 



with 



1 - l \ 1 1 

< a >_ 2 [ ^m + ^} 1 ( ' J>(,) 



Finally using (16), but again with an undetermined constant, k, 



dxG ,.iV c 2 -l 1 1-A . ...o, o 



fc^r<->^^^l (37) 



d 2 M 2 £ ± 4tt 3 A^ c a 2x(A 

as in (31) except for the replacement - 

Thus running coupling effects do not essentially change (30) and (31) 
which now take the forms (34) and (36), respectively. 



6 A simple physical picture 

In this section we present a simple physical picture of p(£±,b, Y) and of 
j^rf- in terms of which gluons in the saturated wavefunction can be viewed 
as determining p and d2 d b x c p i± ■ The picture is complementary to the picture 
given in Ref.9 where d 2^ d 2 b wa s given a simple interpretation in terms of 
scattering of a gluon dipole on a hadron. In Fig. 5 we exhibit the region of 
gluons in saturation using the fixed coupling boundary given in (26). We shall 
carry out the discussion in a fixed coupling context, however, the picture is 
identical for the running coupling case. In Fig. 5 the higher momentum gluons 
are at the top of the figure, with small values of y, while softer gluons are at 
the bottom of the figure with large values of y. In Fig. 6 we plot the x_ extent 
of the gluons in the wavefunction using the Kovchegov gauge where final 
state rescatterings of produced gluons are absent. The x_ boundary of the 
gluons in Fig. 6 is0x_P + w — e y where P + is the momentum of the hadron. 
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We begin our discussion by considering dp/dY which, from (30) and (26), 

is 



dp(£ ± ,b,Y) _iV c 2 -l 2 



dY 7T 



e{. (38) 



The gluons contributing to (38) come directly from parent gluons located in 
the region A of Fig. 5. Region A is at the edge of the saturation boundary and 
is of size A£± ~ £±, Ay ~ 1/ (aN c ). Thus there are on the order of 

(TV 2 — IV 2 

gluons per unit area in A. These gluons emit on the order of ^ N j 2 x • aN c 
gluons per unit area and per unit rapidity directly into the region y ^ Y, k± ^ 
£±. On a scale £± these gluons are essentially uncorrelated in color leading to 



— p{t^Y)~{Nt-l)£l (39) 

agreeing with (38). Two natural questions arise, (i) Why not include emis- 
sions from other areas, such as region B of Fig. 5? Why are there no secondary 
emissions starting from A? That is a source lying in A could emit a gluon at 
transverse momentum on the order of £± and at rapidity y±, lying between 
Y(£±) and Y, and this gluon could then emit a secondary gluon having £± 
and Y. The answer to (i) is that gluons coming from region B are highly 
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shielded on a scale £± and contribute little to charge correlations, although, 
as we shall see shortly they do contribute to the total number of gluons at £± 
and Y. The answer to (ii) is that gluons emitted from yi will typically have 
transverse momentum on the order of Q s (Yi). The probability that a gluon 
be emitted from £±, y x to £+, Y is proportional to £ 2 ± /Q 2 (yi) « 1. Thus the 
gluons that contribute to M^AD are themselves gluons having transverse 
momentum on the order of £± and which come by direct emission from the 
region near £±,Y(£±) where Q S [Y(£±)] = £±. 

Next we turn to (p^it*~d 2 b evaluated at Y. Now gluons coming from sources 
having rapidities between Y(£±) and Y are important. Focus on sources at 
rapidity y x . We can expect a contribution of size 



dxG g|fai) cC-t L 

in — —' a 'QiM im (40) 

coming from source in a range dy\. Of the four factors on the right-hand side 
of (40) the first, Q 2 /a, is the number of gluons per unit area and per unit 
rapidity at y\. The second factor, a, is the strength of emission of gluons 
from y\. d 2 £±/Q 2 is the probability that the emitted gluon end up in the 
momentum range d 2 £±. (Recall that in the Kovchegov light-cone gauge a 
gluon emitted from yi is gauge rotated, locally, by all gluons having larger 
values of a;_H, that is essentially by all gluons having rapidities greater than 
yi as is clear from Fig. 6. These gauge rotations make the typical transverse 
momentum of gluons emitted from sources at y\ on the order of Q s (yi).)dyi 
is the range of sources in rapididty. Integrating dy\ between Y(£±) and Y 
gives!, HI 



^ ^-m)]~^ln^P (41) 



d?£ ± d?b 1 v /J aN c 2 X (\ ) 

which contains the essential factors in (31) and (35). It is important to 
note that the gluons contributing to (39) are highly correlated in color. If 
they were randomly distributed in color they would lead to a contribution 
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X_ = \(2(Xq-X 2 ). 



Figure 6: 



to dp( ^ b ' Y) of size a[Y - Y(£ ± )} ~ ln[Q2(Y)/^J w hi c h is too large. Because 
the sources from which most of the gluons in (41) come from are strongly 
shielded on scale £± the gluons in (41) should themselves contribute little to 
p(£±,b,Y). Finally, gluon sources at y\ having momentum k±_ <C Q s (yi) are 
too few in number to contribute to (40) while gluon sources having kj_ ^> 
Qs{Ui) are too strongly shielded to be important. 

It is interesting to carry out one last calculation which helps explain why 
(28) and (34) have all Y — dependence occurring in the exponent. To that 
end we turn to an interpretation of (25). In the Kovchegov picture the fact 
that is proportional to S comes about because the virtual correction 
shown in Fig. 3 occurs long before the actual scattering, the term shown in 
the figure, or long after the scattering. The exponential result then comes 
about because the virtual gluon corrections are ordered in time, according 
to the longitudinal momentun of the gluons, thus naturally leading to an 
exponential result. How does this occur in a picture where the Y — dependence 
is put into the target rather than in the dipole? When the rapidity of the 
target is increased by dY, there is an increase in the charge correlator by 
dp(£±, b, Y) and this increased charge is located at "large" negative values of 
x_ as shown in the shaded part of Fig. 7 where the top part of that figure 
shows the distribution of charge as a function of Y, much as in Fig. 6. Also 
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Figure 7: 



shown is a scattering off that additional charge by the dipole which reaches 
the charge in dp before it reaches the charge located at higher rapidities. 
There must be at least two gluons exchanged because the charges located at 
different rapidities are too weakly correlated to give a non-zero result and 
there will be no more than two gluons going from dY to the dipole since the 
coupling to the charge in dY is not very large unless dY is quite big. More 
precisely we can rewrite (25) as 



*S, bY) _ f dp(t±,b,Y) d>£ ± g ' 

The factors on the righthand side of (42) are: g 2 /2N c is the coupling of two 
gluons to one of the quarks in the dipole; the 2 counts the number of quarks 
in a dipole while gives the phase space and propagators. S(x±, b, Y) 

then corresponds to the later, in time, scatterings with the rest of the charge 
in the hadron's wavefunction. The integral in (40) goes from £± = l/x± to 
£ ± = Q S (Y). Using (38) then leads to (25). 
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